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EFFECTS  OF  ATOMIC  STRUCTURE  ON  THE  RADIATION 
DYNAMICS  OF  AN  OPTICALLY  THICK  PLASMA 

I.  INTRODUCTION 

The  problem  of  quantitatively  characterizing  the  dynamics  of  high-temperature,  high- 
density  plasmas  is  one  which  has  received  much  attention,  motivated  for  the  most  part,  by  the 
continued  interest  in  inertial  confinement  fusion.  However,  the  complexity  of  the  problems  has 
spawned  many  efforts  which  attempt  to  isolate  one  or  more  of  the  physical  processes  occurring 
in  these  plasmas,  in  particular,  laser-plasma  interaction,  thermal  and  non-thermal  particle  tran- 
sport, plasma  compression,  and  radiation  emission  and  transport,  to  name  a few.  Two  basic 
advantages  of  approaching  the  problem  in  this  manner  are  (i)  the  opportunity  to  analyze  a spe- 
cialized model  of  somewhat  reduced  complexity,  and  (ii)  the  greater  sophistication  used  in 
describing  the  smaller  but  more  involved  subsystem,  both  of  which  ultimately  lead  to  a greater 
understanding  of  the  subsystem  than  would  be  realizable  from  a larger,  more  approximate 
numerical  description  of  the  fully-coupled  radiation-hydrodynamic  problem.  This  report  con- 
cerns itself  with  one  facet  of  the  complex  processes  occurring  in  transiently  heated  dense  plas- 
mas: the  ionization  dynamics  and  radiation  emission,  and  its  effects  on  the  overall  energetics 
and  resultant  diagnostics  associated  with  the  plasma. 

We  consider  a small  volume  element  of  plasma,  homogeneous  in  total  ion  density,  tem- 
perature and  material,  heated  by  uniform  deposition  of  energy,  and  radiating  in  nearly  isotropic 
Manuicript  fubmitted  June  1,  1978. 
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manner.  We  describe,  in  detail,  the  atomic  ionization  state  of  the  plasma,  the  radiation  emis- 
sion resulting  from  atomic  transitions  between  states,  and  the  trapping  of  this  radiation  as  it 
interacts  with  the  plasma.  We  characterize  these  physical  processes  with  a one-dimensional 
model  which  involves  a self-consistent  set  of  time-dependent  equations  describing  the  ioniza- 
tion, radiation,  and  temperature  dynamics  during  the  plasma  evolution.  We  treat  a stationary 
plasma  thereby  isolating  the  ionization  and  radiation  calculation  from  magnetohydrodynamic 
(MHD)  effects  due  to  macroscopic  plasma  motion.  A known  rate  of  energy  absorption  is 
assumed  throughout  the  plasma,  obviating  the  need  for  modeling  an  energy  absorption  mechan- 
ism. 

Since  this  is  not  a coupled  MHD-radiation  calculation,  it  is  our  intent  to  apply  the  model 
to  problems  which  dictate  the  need  for  detail  and  flexibility  in  regard  to  the  physical  processes 
we  have  chosen  to  include.  Although  the  possible  areas  of  study  are  numerous,  we  list  several 
problems  which  are  currently  under  investigation  as  part  of  our  initial  effort: 

1.  Level  structure  — as  a basis  for  accurate  prediction  of  radiation  emission  from  heated 
plasmas,  the  atomic  level  structure  of  the  ions  should  play  an  integral  role.  The  number  of  lev- 
els to  be  included  in  a model,  is,  of  course,  determined  by  the  physics  of  the  problem.  For 
example,  valuable  diagnostic  information  may  be  lost  or  distorted  due  to  neglect  of  some 
atomic  levels,  to  say  nothing  of  the  possible  effect  on  the  gross  energy  balance  of  the  plasma. 
This  model  has  been  designed  specifically  to  allow  total  flexibility  in  the  selection  of  the  atomic 
structure.  For  purposes  of  comparison,  ionization  stages,  energy  levels,  and  :adiation  transport 
may  be  easily  and  quickly  altered  to  facilitate  our  structure  investigations.  One  aspect  of  this 
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work  will  be  to  determine  the  error  incurred  in  the  radiation  emission  due  to  the  neglect  of  cer- 
tain atomic  levels  in  various  materials,  thereby  establishing  standards  for  completeness  with  a 
minimum  essential  amount  of  structural  detail,  to  be  used  in  further  investigations. 


2.  Isoelectronic  Sequences  — an  important  part  of  the  level  structure  problem  is  the  predic- 
tion of  transition  rates,  energy  levels,  and  radiation  emission  (as  a function  of  a given  energy 
level  or  groups  of  levels)  in  similarly  charged  ions  as  a function  of  nuclear  charge.  These 
trends  should  be  used  to  establish  different  requirements  for  adequate  ionization  — radiation 
descriptions  as  one  goes  from  moderate  to  high-Z  materials. 

3.  Rale  Coefficients  — the  validity  of  an  atomic  level-structure  model  depends  fundamen- 
tally on  how  well  the  transitions  between  these  levels  are  described.  It  is  common  practice  to 
employ  general  prescriptions  to  calculate  rate  coefficients  that  are  independent  of  the  individual 
characteristics  of  a specific  ionic  transition.  While  having  the  advantage  of  simplistic  and 
economic  implementation,  the  accuracy  of  these  coefficients  may  be  insufficient  for  the 
intended  model.  The  alternative  is  to  treat  each  transition  separately  and  calculate  accurately 
(but  at  greater  expense)  the  rate  coefficient  over  a temperature  range  of  interest.  A problem 
presently  under  investigation  is  to  assess  the  errors  incurred  when  approximate  rate  coefficients 
are  introduced  into  calculations. 

4.  Population  Inversions  — an  interesting  and  useful  application  of  a detailed  atomic  model 
is  in  the  area  of  X-ray  and  VUV  plasma  dynamic  lasers.  The  calculation  can  be  used  to  investi- 
gate energy  level  pairs  which  appear  promising  as  inversions  when  an  ionized  plasma  is  allowed 
to  expand,  cool  and  recombine  Wc  expect  the  knowledge  gained  from  the  study  of  isoelec- 
tronic sequence  trends  in  level  populations  will  play  a significant  part  in  predicting  how  the  gain 
from  an  inversion  will  scale  with  atomic  number. 
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S.  Diagnostics  — although  a given  model  may  neglect  some  effects  for  the  sake  of  simpli- 
city and  place  specific  emphasis  on  a few  processes,  it  should  nevertheless  be  capable  of  provid- 
ing predictions  for  experiment.  For  example,  the  calculated  line  emission  spectra  can  be  com- 
pared to  experimentally  obtained  data  as  a measure  of  the  temperature  and  density  in  the 
plasma  as  a function  of  heating  conditions  and  plasma  atomic  number. 

II.  MODEL  DESCRIPTION 


The  theoretical  description  of  our  model  is  logically  divided  into  three  parts:  0)  the 
plasma  heating,  thermal  conduction,  and  temperature  calculation,  (2)  the  atomic  structure  and 
ionization  dynamics,  and  (3)  the  radiation  transport  model.  The  first  and  third  topics  are  com- 
mon to  all  calculations,  while  the  second  is  particular  to  the  specific  plasma  under  investigation. 


A.  Energy  Deposition  and  Transfer  — the  details  of  the  basic  plasma  physics  included  in 
our  model  are  described  elsewhere1;  only  the  basic  equations  are  presented  here.  The  set  of 
equation  solved  include  an  electron  and  ion  temperature  equation. 


< 


l 


_a_ 

bt 


y N'  kB  T, 


a. 

bt 


bE,  «. 

+ — - S'  - 
bt 


3«,  Tt  _ 3 A/,  N 'J± 
y I,  M,  r,  #' 


f n,  kB  r,  - 


3A/,  N'k„ 

M,  r. 


IT,-Tt 


e) 


- R 


(1) 

(2) 


where  A/,.  A/,,  T„  T,  and  Nt,  N,  are  electron  and  ion  mass,  temperature  and  density  respec- 
tively, k,  is  electron  thermal  conductivity1,  kB  is  Boltzmann’s  constant,  5,  is  the  energy  source 
term,  r,  is  the  electron-ion  collisional  equipartition  time1,  A is  a term  characterizing  radiation 
loss  by  the  plasma,  and  £,  is  energy  expended  by  the  electrons  during  ionization  or  excitation 
processes.  The  term  /,  represents  a temperature  gradient  scale  length  and  is  determined  in  such 
a way  as  to  realistically  model  conduction  losses  to  regions  of  the  plasma  external  to  the  sphere. 
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The  radiation  loss  term  R , contains  contributions  corresponding  to  bremsstrahlung,  radiative 
recombination,  and  line  emission  processes  and  will  be  defined  more  fully  in  the  next  section. 
The  source  term  S,.  can  model  energy  deposition  at  any  desired  rate  or  pulse  shape;  presently  a 
Gaussian  rate  is  being  used. 


B.  Atomic  Physics  — the  method  of  solution  of  the  rate  equations  is  identical  for  any 
atomic  system;  only  the  number  of  equations  to  be  solved  will  change.  The  model  solves  a set 
of  time-dependent  rate  equations  containing  terms  which  describe  the  population  and  de- 


population of  atomic  levels.  The  equations  are  integrated  in  time  using  an  iterative  predictor- 
corrector  algorithm  to  obtain  convergence  while  avoiding  "stiffness"  by  separate  treatment  of 
asymptotic  rates.  The  method  of  calculation  of  the  rate  coefficients  is,  of  course,  fundamental 
to  the  accuracy  of  our  ionization  and  radiation  calculation  and  it  is  worthwhile  to  discuss  several 
methods  in  detail.  The  processes  presently  included  ir.  our  model  are  collisional  ionization  and 
recombination,  radiative  recombination,  collisional  excitation  and  de-excitation,  spontaneous 


emission,  and  stimulated  emission  and  absorption.  The  first  three  processes  have  been  dis- 
cussed elsewhere1,  while  the  last  two  will  be  discussed  in  the  next  section;  the  spontaneous 
emission  is  given  by  the  Einstein  A coefficient4.  The  collisional  de-excitation  process  is 
obtained  from  detail  balancing  of  the  collisional  excitation  rate;  the  latter  has  been  calculated  in 
four  ways.  For  approximate  collisional  rates,  a modified  version  of  Van  Regemorter's  semiclas- 
sical  method5  has  been  used. 
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where  lH  is  the  Rydberg  energy,  A£  is  the  transition  energy,  /„  is  the  transition  oscillatoi 
strength,  and  g,,  is  an  average  Gaunt  factor.  In  this  approximation,  Allen’s  values*  have  beer 
used  and  parameterized  by7 
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<g>  - A + (BY  + C)  In 


where  A,  B , and  C are  constants,  and  Y - \E/Tr. 
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When  greater  accuracy  is  required,  the  rates  are  calculated  by  the  method  of  distorted  waves.1 
We  have  calculated  rates  for  several  hydrogenlike  and  heliumlike  ions  along  the  isoelectronic 
sequence  using  the  Born  approximation  as  parameterized  by  Vainshtein  and  Sobel’man9, 
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where 


Here  E is  the  ionization  energy,  / is  the  orbital  angular  momentum,  0 - A£/kfl  Tr,  and  a and 
X are  parameters  characteristic  of  the  transition.  Lastly,  in  the  case  of  Al  XIII  (hydrogenlike) 
we  have  taken  several  transition  rates  from  Landshoff  and  Perez10,  which  were  calculated  using 
a Coulomb-Born  approximation  (without  exchange).  A comparison  of  these  methods  is  made 
in  a following  section. 


The  atomic  structure  has  been  developed,  at  present,  for  three  plasmas,  carbon,  alumi- 
num, and  argon.  The  carbon  and  aluminum  structure  is  shown  in  Fig.  1 while  the  argon  model 
is  depicted  in  Fig.  2.  Excited  states  are  coupled  to  neighboring  ground  states  by  collisional  ioni- 
zation and  recombination  (above)  and  collisional  and  stimulated  excitation  and  de-excitation 
and  spontaneous  emission  where  allowed  (below).  In  addition,  the  excited  states  are  coupled  to 
each  other  by  the  last  $ processes  for  all  ions  but  Ar  XII  and  Ar  XIV.  All  adjacent  ground 


states  are  coupled  by  the  first  3 processes. 
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C.  Radiation  Transport  — a phenomenological  radiation  transport  scheme  has  been  added 
to  the  existing  formalism  to  extend  the  predictive  capability  of  the  calculations  to  cases  where 
the  plasma  is  optically  thick.  The  model  used  is  fundamentally  one  developed  by  Julienne  and 
Davis"  and  is  based  on  photon  escape  from  the  plasma  by  frequency  diffusion  into  the  wings  of 
the  broadened  line  profile.  Effective  rates  for  photon  scattering  by  ions  and  photon  escape  rates 
have  been  calculated  and  are  given  by. 
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where  t„  is  the  optical  depth  of  the  plasma'*  at  line  center. 
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X)  is  the  frequency  (in  Doppler  widths)  at  the  point  in  the  line  profile  where  the  optical  depth 
drops  to  1,  L is  the  effective  path  length,  c is  the  speed  of  light,  g is  the  statistical  weight,  and 
N,  and  N,  are  the  population  densities  of  the  lower  and  upper  states,  respectively.  The  Doppler 
width  is  given  by. 
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where  tr0  is  the  line  center  frequency  and  M is  the  ion  mass.  For  a Doppler  broadened  line 
profile,  x j - (In  t„),/j,  thus,  giving  r,  - 2 In  r„  L/c  and  r,  - 2(ln  t()i/!  L/y/n  t„c.  We  have 
also  extended  the  calculation  to  include  Voigt  line  profiles  by  using  .a  approximate  expression 
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where,  the  damping  parameter  a is  given  by. 
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Accurate  Voigt  and  Doppler  profiles  have  been  computed  numerically  to  test  the  validity  of  the 
above  expression,  and  the  results  are  shown  for  several  lines  of  aluminum  K -shell  transitions  in 
Figs.  3-5.  A complete  tabulated  set  of  results  for  various  temperatures  and  densities  is  given  in 
Table  I,  where  x,  (theory)  is  set  by  the  maximum  of  the  Doppler  or  Voigt  values  as  given 
above  and  compared  with  an  accurately  calculated  X|  (numerical)  taken  from  the  actual  profiles 

(L  — 0.01  cm). 


Separate  rate  equations  are  now  set-up  for  each  transition  describing  the  density  ot  pho- 


tons in  the  plasma. 
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where  j and  / denote  upper  and  lower  states  and  A is  the  Einstein  coefficient.  The  terms  on  the 
right  hand  side  characterize,  in  order,  the  processes  of  spontaneous  emission,  stimulated  emis- 
sion, stimulated  absorption,  and  escape  from  the  plasma.  Each  rate  equation  for  state  densities 
will  also  be  modified  by  the  transport  terms  (where  appropriate  transitions  occur)  giving, 

“ I w„  N,  - £ W„  N,  + j-N;  Nr  <U) 

where  W represents  the  processes  referred  to  in  the  previous  section. 
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In  practice  equations  5 and  6 become  undefined  at  small  values  of  r,„  hence,  smooth 
interpolations  were  added  to  keep  the  transition  from  an  optically  thick  to  an  optically  thin  line 
analytically  defined. 

This  radiation  transport  model  relies  on  the  fact  that  the  photon  distribution  from  a given 
transiton  will  spread  sufficiently  into  the  wings  to  allow  stabilization  on  a time  scale  short  com- 
pared with  plasma  time  scales  of  interest.  This  equilibration  time  is  given  by" 
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Since  typical  plasma  lengths  are  less  than  10~2  cm  and  time  scales  for  the  rate  equations  are 
~ 10‘12  sec,  f|  is  relatively  smaller  for  even  marginally  thick  lines.  That  is  to  say,  as  a line 
becomes  thick  enough  to  require  a transport  treatment,  the  line  profile  has  already  broadened 
sufficiently  for  equilibration  to  satisfy  the  above  criteria. 

The  radiation  calculation  to  determine  the  plasma  emission  is  accomplished  using  three 
separate  equations,  one  each  for  bremsstrahlung12  (with  unit  Gaunt  factor), 
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where  a is  the  radiative  recombination  rate  coefficient.  The  three  contributions  are  then 
summed  and  inserted  in  the  electron  temperature  equation  as  the  term  R. 

III.  PRELIMINARY  INVESTIGATIONS 

In  preparation  for  the  studies  discussed  in  section  II,  we  have  used  the  code  to  shed  light 
on  some  interesting  preliminary  problems  as  well  as  to  compare  our  results  with  other  similar 
but  more  detailed  calculations.  Some  of  these  results  will  be  presented  in  the  following  sec- 
tions. 


A.  Comparison  of  Rale  Coefficients  — several  methods  have  been  used  to  calculate 
electron-ion  collisiona!  excitation  rate  coefficients  as  was  mentioned  in  section  II.  Comparisons 
have  been  made  between  these  methods  for  the  purpose  of  weighing  the  inaccuracy  introduced 
into  the  simulation  by  using  more  approximate  (but  relatively  simpler)  calculations  for  the 
rates.  The  first  comparison  studied  the  results  of  a distorted  wave  (DW)  calculation,  a semi- 
classical  impact  parameter  (SCI)  approximation,  and  a parameterized  Born  approximation 
(BORN)  as  applied  to  determine  several  collisional  excitation  rate  coefficients  of  the  heliumlike 
ion  of  aluminum.  The  BORN  method  was  used  to  calculate  rate  coefficients  from  two  different 
sets  of  parameters,  one  describing  hydrogenic  systems'0  and  one  for  heliumlike  species" 
(designated  BORN(//>  and  BORN(We)).  The  results  are  depicted  graphically  in  Figs  6-14.  The 
SCI  method  is  shown  only  for  allowed  lines,  since  no  simple  semi-classical  prescription  is  avail- 
able for  forbidden  transitions.  In  addition,  the  BORN  method  is  shown  only  for  those  transi- 
tions which  were  specifically  parameterized  by  the  authors'"  " (hence,  no  An  — 0 transitions 
for  BORN(W)  and  only  n-1  ton-2andn-2.  AN  - 0 transitions  for  BORN  (Me)).  Since 
we  consider  the  DW  calculation  to  be  the  most  sophisticated  of  those  considered  here,  we  will 
use  it  as  the  standard  of  accuracy  by  which  we  make  comparisons. 
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Consider,  first,  the  SCI  and  DW  calculations.  For  Aa  « 1 transitions,  the  quantitative 
agreement  appears  som.  what  good,  however,  qualitatively,  the  temperature  dependence  of  the 
SCI  rate  coefficients  is  clearly  incorrect  at  higher  temperatures  The  accuracy  of  the  method 
deteriorates  as  the  transition  energy  decreases;  notice  the  poor  agreement  between  the  3 d - 4 p 
and  2 p - 3 d transitions  (Fig  9)  and  the  Aa  — 0 transitions  (Fig  12).  The  agreement  is  rea- 
sonably accurate  for  transitions  from  the  ground  state  for  which  A/  - 1.  The  accuracy  of  the 
SCI  method  is  strongly  dependent  on  the  value  of  the  oscillator  strengths  employed  in  the  cal- 
culation. To  obtain  oscillator  strengths  for  higher  Z materials,  we  extrapolated  many  tabulated 
values  from  lower  Z ions4  and  compared  the  results,  quite  favorably,  with  studies  done  by  Lin 
et  at. 14  Schiff  et  at, 15  Weiss,"1  and  Smith  and  Wiese;11  the  results  are  shown  for  a few  transi- 
tions in  Figs.  15-18. 

Definite  trends  are  evident  when  the  BORN(//)  method  is  compared  to  the  DW  calcula- 
tion quantitatively.  For  large  energy  separations  the  BORN(//)  formula  yields  consistently 
lower  values  for  the  excitation  rate  by  about  30%,  as  typified  by  the  the  Lyman  transitions. 
This  trend  reverses  itself  for  the  smaller  energy  differences  such  as  the  Balmer  and  Paschen 
series,  where  the  BORN(//)  calculation  gives  typically  greater  values.  The  parameter  constants, 
a and  \ used  in  formula  4 were  calculated  for  a hydrogenic  system,  however,  and  the  agree- 
ment between  these  two  methods  is  expected  to  improve  with  increasing  Z as  the  heliumlike 
ion  more  closely  resembles  a single  electron  orbiting  a highly-charged  core. 

Several  rate  coefficients  were  calculated  using  the  BORN  approximation  but  employing 
parameters  fitted  to  heliumlike  transitions;1’  the  results  ‘ave  been  plotted  in  Figs.  6,7,12,13,14 
for  a few  transitions,  including  some  involving  change  of  spin.  The  Is  - 2s1  transition  com- 
pares excellently  with  the  DW  rate,  and  the  Is  - 2 />’  also  gives  a favorably  representation,  but 
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the  temperature  dependence  of  the  Is  — 2/» 1 transition  is  not  accurate  at  higher  temperatures 
(Fig.  7).  In  addition,  the  agreement  between  the  two  methods  is  very  poor  for  the  2s1  - 2 p\ 
2 p’  - 2s',  and  2 p}  - 2 p'  transitions  (Figs.  12  and  14),  both  quantitatively  and,  as  seen  in  Fig. 
14,  qualitatively. 

A second  rate  coefficient  comparison  was  made  between  the  BORN(W)  and  DW  methods 
when  applied  to  the  hydrogenlike  ion;  also,  the  rates  calculated  independently  using  a 
Coulomb-Born  approximation  without  exchange10  were  compared  where  possible  to  check 
agreement  with  our  hydrogenlike  aluminum  rate  coefficients.  The  DW  and  BORN  calculations 
were  checked  for  C VI  as  well  as  Al  XIII  to  establish  whether  the  accuracy  of  the  BORN 
method  was  Z-dependent.  In  all  cases  calculated,  the  discrepancy  between  the  three  methods 
never  exceeded  25%  and  was  usually  on  the  order  of  5-10%,  with  slightly  better  (but  not 
significant)  accuracy  obtained  for  the  Al  XIII  ion. 

B.  Validity  of  a Statistically  Populated  Hydrogenlike  Ion  — in  order  to  calculate  level  popula- 
tions and  radiation  more  economically,  we  have  represented  the  hydrogenic  ions  in  our  model 
as  a system  of  V levels,  as  was  stated  in  the  previous  section,  rather  than  treat  each  nl  level 
separately.  This  approximation  assumes  that  the  collisional  coupling  between  different  I sublev- 
els of  a given  n level  is  sufficiently  rapid  that  the  sublevels  relax  to  their  statistical  distributions 
on  an  instantaneous  time  scale  relative  to  the  other  transition  times  populating  the  sublevels. 


\ 


In  order  to  verify  this  assumption  we  used  the  condition  derived  by  Sampson1*  to  deter- 
mine the  lowest  n value  at  which  the  states  should  be  statistically  populated, 
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1.1  x >0  '*  — N,  > I (11  ^ 7 < 19)  (18) 

Z" 

with  T,.  in  degrees  kelvin.  For  u carbon  plasma  at  300  eV,  an  electron  density  of  only 
4.3  x 1017  cm  3 is  required  for  equilibration  of  the  n - 2 levels,  a density  which  is  usually 
exceeded  for  plasmas  of  interest  in  this  work.  For  aluminum  and  argon,  the  criterion  yields 
1 x I0-’1  and  1.8  x I0:J  cm  ’ for  n — 2,  and  3.4  x 10’“  and  5.6  x 10 70  cm  1 for  n - 3,  respec- 
tively. The  condition  for  statistical  equilibrium  is  met  easily  for  the  n - 3 levels  by  both  plas- 
mas, since  our  study  is  concerned  with  plasmas  of  ion  density  5 x 10N  cm  3 and  greater  How- 
ever, the  n — 2 levels  are  marginally  unequilihrated , according  to  condition  (18),  and  should 
actually  be  treated  as  separate  levels  The  criterion  may  be  relaxed  somewhat  due  to  the  fact 
that  it  is  derived  for  a plasma  which  is  optically  thin,  whereas  in  reality,  radiation  trapping  due 
to  opacity  in  plasmas  at  these  densities  should  drive  the  level  populations  closer  to  their  LTE 
values.  In  particular,  the  (2.v,  2 />)  level  system  of  the  hydrogenlike  ion  may  not  be  statistically 
equilibrated  due  to  the  depletion  by  strong  radiative  decay  of  the  2 />  state.  In  an  optically  thick 
plasma,  however,  photo-excitation  will  exceed  collisional  excitation  as  the  dominant  mechanism 
populating  this  state,  preferentially  increasing  its  density  compared  to  the  2s  level  and  restoring, 
to  a degree,  the  statistical  equilibrium  of  the  n - 2 level.  Since  the  L„  line  is  consistently  opti- 
cally thick  for  the  plasma  densities  considered  in  this  study,  we  have  chosen  to  represent  the 
n — 2 level  as  a single  degenerate  state  as  we  have  done  with  the  n > 3 levels. 

C.  Degeneracy  of  levels  of  the  helinmlike  ion  — the  structure  of  the  heliumlike  ion  differs 
from  that  of  the  hydrogenlike  ion  in  that  (i)  two  distinct  spin  systems  (ortho-  and  parahelium) 
are  present,  and  (ii)  the  energy  differences  between  the  / components  of  an  n level  are  larger, 
particularly  at  low  n values.  This  would  seem  to  indicate  that  a similar  attempt  to  model  only 
"degenerate"  n levels  of  the  heliumlike  ion  by  assuming  that  the  various  / components  of  both 
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spin  systems  equilibrate  instantaneously  to  their  statistical  populations  would  yield  incorrect 
results.  However,  as  principal  quantum  number  increases,  the  energy  separation  between  /sub- 
levels  of  each  spin  system  may  decrease  sufficiently  that,  once  again,  ion-ion  collisions  may  be 
the  dominant  transition  mechanism  causing  population  changes  among  these  levels.  Moreover, 
it  is  well  known  that  resonance  transitions  occurring  in  the  singlet  system  of  the  heliumlike  ion 
tend  to  be  optically  thick  for  dense  plasmas,  due  to  the  persistently  high  concentration  of  the 
ground  state  of  this  closed-shell  configuration  over  a wide  temperature  tange  Hence,  as  was 
stated  in  the  previous  section,  opacity  effects  should  optically  pump  some  of  the  lower  depleted 
levels,  restoring  their  populations  and  driving  them  closer  to  their  statistical  distributions. 


In  order  to  quantitatively  assess  these  effects  on  the  excited  level  population:  of  the  heli- 
umlike ion  and  to  detemine  at  what  value  of  n,  the  / levels  relax  to  their  statistical  populations, 
we  calculated  the  ion-ion  collision  rates  coupling  the  nl  components  for  the  n - 2.3  levels  of 
the  argon  heliumlike  ion  and  compared  them  with  the  fastest  radiative  rates  depleting  their 
population  densities.  To  calculate  the  rate  coefficient  for  collisional  excitation  by  ions  we  used  a 
formula  by  Drawin19  with  a guillotine  factor  given  by  Weisheit20, 
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(19) 


where  m is  the  reduced  mass  of  the  incident  and  target  ion  and  Z is  the  charge  of  the  incident 
ion.  For  the  n - 2 level,  the  fastest  radiative  rate  was  for  the  (Is  - 2p')  transition, 
1.07  x 10u  sec  while  the  combined  electron  and  ion  excitation  rate  coupling  the  2s1  to  the 
2p'  state  was  2.3  x 10"  sec  '•  at  Tf  - 1.5  keV  and  Nr  - 10JI  cm  \ Clearly,  the  2/>'  level  will 
be  depleted  so  rapidly  by  spontaneous  decay  that  equilibration  with  (he  2s1  will  never  occur  at 
these  plasma  parameters.  It  is  also  very  unlikely  that  opacity  effects  will  repopulate  the  2 p' 
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level  sufficiently  to  overcome  the  3-orders-of-magnitude  difference  between  these  two  rates. 
The  2s]  - 2 p3  collisional  excitation  rate  is  basically  identical  to  the  singlet  transition  (the 
electron-ion  rate  is  larger  but  the  ion-ion  rate  is  smaller),  but  the  radiative  decay  rate  of  the  2 p3 
state  is  only  6.07  x 10"  sec'1.  For  the  assumption  of  statistical  equilibration  to  be  valid,  one 
should  expect  that  the  collisional  rate  should  exceed  the  radiative  rate  by  at  least  a factor  of  ten. 
Although  the  rates  are  close,  opacity  effects  should  be  small  for  the  intercombination  line,  due 
to  the  small  Einstein  A value,  hence  these  levels  should  not  equilibrate  at  the  assumed  plasma 
parameters  either. 


Repeating  the  calculation  for  the  n - 3 level,  we  find  that  the  collision  rate  for  the 
(3s1  - 3d')  transition  will  equal  the  fastest  radiative  decay  rate  (2 p'  - 3d')  at 
Ne  - 7.7  x 1020  cm'5  and  the  same  condition  occurs  for  the  (.3d'  - 3 p')  collisional  rate  and 
the  (Is  - 3 p')  radiative  decay  at  Ne  - 6 x 10JO  cm'3  (level  energies  were  taken  from  Ermolaev 
and  Jones").  At  Nf  ==  10"  cm'3,  conditions  for  statistical  equilibration  of  the  n - 3 level  are 
marginally  satisfied,  but  even  small  opacity  effects  should  have  an  influence  on  the  repopulation 
of  radiatively  depleted  / levels.  An  approximate  measure  of  this  effect  is  the  ratio  of  the 
photo-excitation  to  the  collisional  excitation  rate  for  one  of  these  states.  We  can  express  this 
ratio,  based  on  our  frequency  diffusion  transport  model,  as 


*2.  _ yd  ~ 
x,i  1 + pqy 


(20) 


where  y - r„  (ir  In  r„),/J  and  p9  is  the  quenching  parameter  defined  by  p9  — Y/(A  + K),  with 
A being  the  Einstein  A of  the  line  and  Y the  rate  accounting  for  all  other  processess  which 
deplete  the  level,  excepting  A.  Typical  values  of  these  quantities  at  the  plasma  conditions  pre- 
viously quoted  are  p„  — 0.3  and  r0  — 15  (for  a 100  micron  plasma),  yielding  a ratio  of  2.17  for 
the  (Is1  — 3/»’>  transition.  Hence,  doubling  the  rate  at  which  the  ground  state  populates  the 
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ip'  state  should  have  a significantly  beneficial  effect  on  equilibrating  the  n - 3 singlet  levels. 
The  triplet  system  for  this  level  should  parallel  the  singlet  system  for  the  same  reasons  which 
were  stated  in  the  discussion  of  the  n - 2 system. 

In  view  of  the  above  study,  we  have  chosen  to  include  the  four  individual  nl  components 
in  our  model  for  the  n - 2 level  system,  hence,  taking  explicit  account  of  all  processes  coupling 
each  of  the  four  states.  However  based  on  our  discussion  of  the  n - 3 level  we  expect  that 
collisions  will  dominate  the  higher  states,  sufficiently  aided  by  pumping  due  to  radiation  trap- 
ping, and  we  have  chosen  to  describe  all  higher  levels  (n  3)  of  the  argon  heliumlike  ion  as 
two  degenerate  n states,  a singlet  and  a triplet.  In  the  argon  model  used  in  our  preliminary 
investigations,  however,  we  have  omitted  the  n -3. 4.  and  5 triplet  states  as  an  initial  approxi- 
mation. 

D.  Comparison  with  a Similar  Model  — We  have  compared  the  results  predicted  by  our 
model  with  that  of  a similar  calculation22  designed  for  implementation  in  cylindrical  geometry 
calculations.  However,  by  analyzing  calculations  of  the  atomic  state  of  the  plasma  in  the  opti- 
cally thin  approximation,  meaningful  comparisons  can  be  made  between  the  two  models.  Both 
simulations  were  at  constant  ion  density  of  5 x lO1’  cm  } and  electron  temperature  of  550  eV 
for  an  aluminum  plasma.  The  only  differences  in  the  two  calculations  were  in  the  structure  of 
the  K-shell  excited  levels,  and  some  relevant  conclusions  have  resulted  from  this  study. 

The  heliumlike  level  system  was  nearly  identical  in  both  models,  with  the  exception  that 
the  5 p]  level  was  missing  in  our  calculation  (see  Fig.  1).  The  populations  of  the  included  states 
are  shown  in  Table  11  for  both  simulations.  The  differences  between  the  densities  are  less  than 
15%  in  all  cases  except  the  Is2/>J  level  where  the  error  is  37%.  This  is  attributable  to  the  fact 
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that  slightly  different  decay  rates  were  used  for  the  (Is  - 2 />')  transition”  in  the  two  models. 
The  consistency  of  both  calculations,  however,  is  clearly  evident. 

Contrasting  differences  are  noticeable  in  the  results  obtained  for  the  hydrogenlike  ion. 
Our  model  contained  a degenerate  system  of  the  lowest  3/r  levels,  with  each  / component 
assumed  to  be  populated  statistically  (see  last  section),  while  the  other  calculation  included  only 
p states  up  to  the  4 />'  level.  In  addition,  no  collisional  coupling  processes  were  included 
between  these  p levels.  The  densities  shown  in  the  table  for  p states  from  our  model  are 
obtained  from  the  n-level  densities  by  simple  degeneracy  ratios.  Significant  differences  are  seen 
to  occur  between  the  two  models,  becoming  larger  as  the  principal  quantum  number  increases. 
We  attribute  this  trend  in  the  differences  to  the  progressively  smaller  fraction  of  the  total  n 
population  occupying  the  p state  as  n increases  above  2,  indicating  the  necessity  for  including  a 
more  comprehensive  level  structure  for  the  lower  n values. 

E.  Comparison  of  Optically  Thick  and  Thin  Calculations  — To  see  the  profound  effect  of 
radiation  trapping  on  the  calculation,  two  separate  runs  modeling  a laser-heated  plasma  were 
made  with  identical  input  parameters,  changing  only  the  transport  calculation  from  perfectly 
thin  to  thick.  The  results  presented  are  for  a 77  micron  radius  sphere  of  aluminum  plasma  at  a 
density  of  7.6  x 10”  ions/cm3.  The  laser  pulse  was  0.2  nsec  (FWHM)  and  a total  of  0.8  Joules 
of  laser  energy  were  absorbed  in  the  sphere.  Electron  temperatures  peaked  along  with  the  laser 
pulse  at  about  430  eV  and  decayed  slowly  to  about  110  eV.  The  line  spectra  for  the  thick  case 
are  shown  in  Figs.  19  and  20,  with  solid  horizontal  bars  indicating  the  emission  in  the  thin  case. 
In  general,  the  optically  thick  lines  will  decrease  and  thin  ones  will  increase  when  radiation  is 
trapped.  The  L„  and  lines  of  Al  XU  are  seen  to  decrease  since  they  are  thick,  while  all 
other  Al  XII  lines  increase  due  to  the  enhanced  population  from  the  upper  levels  of  the  thick 
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transitions  causing  greater  overall  excited  state  populations.  However,  all  or  the  Al  XIII  lines 
are  seen  to  increase  over  the  thin  case,  even  though  the  L„  line  is  marginally  reabsorbed 
(r „{MAX)  ~ 6).  The  apparent  increase  is  due  predominantly  to  the  pumping  of  the  Al  XIII 
ground  state  by  collisional  ioni/ation  of  the  Al  XII  excited  states,  now  with  increased  popula- 
tions. The  Al  XIII  ground  state  now  couples  to  its  excited  states  in  turn,  increasing  their  popu- 
lations. To  verify  this,  another  run  was  made,  ("thick-thin")  with  the  heliumlike  lines  tran- 
sported but  the  hydrogenlike  lines  untransported.  The  results  are  shown  on  the  graphs  by  the 
dotted  horizontal  lines  on  the  Al  XIII  lines  (Al  XII  line  emission  was,  as  expected,  unchanged). 
Note  the  increase  in  the  H(L„ ) line  in  the  "thick-thin"  case  over  the  thin  case  and  then  another 
decrease  as  both  ions  are  calculated  as  thick;  the  result,  therefore,  is  a net  increase  in  the  emis- 
sion of  this  line  even  though  it  is  somewhat  thick 

The  effect  of  photon-pumping  on  the  upper  levels  of  thick  lines  has  a profound  effect  on 
the  thin-line  emission;  order  of  magnitude  increases  are  seen  to  occur  in  H(L„).  H,„  //„.  He,, 
and  Hefi  line  emission  corresponding  to  similar  increases  in  the  upper  states  of  these  transi- 
tions. The  heliumlike  intercombination  line  also  displays  some  enhancement  since  the  2 />'  state 
is  strongly  collisionally  coupled  to  the  pumped  2/>'  state  at  these  electron  densities.  It  is  hoped 
that  some  of  the  more  dramatic  "thickness"  effects  discussed  here  might  be  exploited  for  diag- 
nostic purposes  to  determine  spatially  resolved  parameters  in  these  plasmas. 

IV.  DISCUSSION 

In  an  attempt  to  shed  some  light  on  a few  of  the  more  subtle  problems  that  exist  in  the 
mathematical  modeling  of  the  atomic  state  of  high-temperature,  high-density  plasmas,  we  have 
been  able  to  reach  some  significant  conclusions  based  on  our  investigations.  In  addition,  the 
results  obtained  thus  far  will  guide  our  future  attempts  to  model  atomic  structure  and  its  radia- 
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tive  effects  on  the  energy  balance  of  these  plasmas,  and  therefore,  in  summary,  a disussion  of 
these  results  is  fruitful. 

The  acquisition  of  reliable  rate  coefficients,  in  particular,  those  describing  collisional  exci- 
tation and  de-excitation,  have  always  presented  a major  obstacle  to  those  who  would  include 
ionization  dynamics  in  their  plasma  models.  Many  of  the  scattering  methods  are  too  detailed 
and  involve  expensive  computer  codes  to  generate  rates  which  are  calculated  with  more  accu- 
racy than  the  plasma  model  may  initially  demand.  Other  simpler  calculations  may  not  give  a 
comprehensive  treatment  of  all  rates  of  a given  ionization  stage,  or  are  based  on  assumptions 
which  introduce  intolerable  or,  more  often,  unknown  inaccuracies  into  the  rates  In  view  of 
our  preliminary  results  in  this  area,  we  feel  several  conclusive  statements  can  be  made  The 
well-known  semi-classical  impact  parameter  formula  can  be  used  to  generate  rates  for  dipole 
allowed  transitions  with  varying  degrees  of  error.  For  most  resonance  transitions  of  the  heli- 
umlike system  the  oscillator  strengths  have  already  saturated  at  medium  / to  near  their  hydro- 
genic  values,  and  can  be  obtained  by  extrapolation  techniques  to  good  accuracy,  thereby  yield- 
ing rates  accurate  to  about  10%.  with  our  standard  of  comparison  being  set  by  a reliable 
distorted-wave-with-exchange  calculation.  For  transitions  of  lower  energy  difference,  this 
method  becomes  unreliable,  yielding  large  quantitative  differences  along  with  a more  deviant 
temperature  dependence.  The  values  generated  for  transitions  in  which  An  - 0 may  he  as 
much  as  an  order  of  magnitude  in  error,  making  this  a totally  unacceptable  method  for  this 
case.  In  addition,  the  method  is  undefined  for  forbidden  transitions.  A simple  hydrogenic 
Born  approximation,  on  the  other  hand,  yields  rates  for  all  transitions  within  a heliumlike  spin 
system,  and  they  tend  to  agree  favorably  with  the  temperature  dependence  of  the  more  accurate 
DW  rates.  The  error  is  larger  for  resonance  transitions  than  it  is  with  the  SCI  calculation  but 
yields  better  agreement  for  other  higher  level  transitions. 
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Specific  parameters  developed  for  heliumlike  ions  were  also  calculated  in  the  Born  approx- 
imation and  found  to  give  very  inconsistent  results.  Rates  for  some  transitions  differed  by  less 
than  10%  from  the  DW  calculation  while  others  showed  discrepancies  of  an  order  of  magnitude 
or  more  with  a totally  dissimilar  temperature  dependence  from  that  predicted  by  the  distorted 
wave  calculation. 

However,  encouraged  by  the  excellent  agreement  between  the  BORN  and  DW  calcula- 
tions for  several  transitions,  it  is  conceivable  that  a complete  set  of  collisional  rates  could  be 
calculated  for  heliumlike  ions  using  the  BORN  formula,  equation  (18),  if  a more  appropriate 
and  comprehensive  set  of  a and  x parameters  could  be  determined,  and  that  these  parameter 
fits  would  compare  quite  favorably  with  the  rates  calculated  using  more  accurate  but  more 
detailed  scattering  theories. 

Although  we  have  sought  to  investigate  simple  methods  of  generating  an  approximate  set 
of  rate  coefficients  for  use  in  a general  K-shell  structural  model,  we  do  not  wish  to  underem- 
phasize the  value  or  accuracy  of  a distorted  wave  calculation,  specifically  for  An  — 0 and  spin- 
flip  transitions,  nor  the  importance  of  its  use  in  benchmarking  any  attempt  at  an  empirically 
determined  parameter-fit  of  rate  coefficients  as  we  have  done  here. 

In  the  case  of  the  hydrogenlike  ion,  we  have  established  that  the  Born  approximation 
yields  excellent  agreement  with  a DW  treatment  of  collisional  excitation  rates,  and  that  the 
agreement  improves  with  increasing  Z,  yielding  a simple  and  economical  method  of  obtaining  a 
complete  set  of  accurate  values.  In  addition,  we  have  shown  that  under  certain  plasma  condi- 
tions, little  or  no  accuracy  is  lost  in  employing  a single  level  representation  for  the  various  / 
states  of  a given  n level  for  this  ion,  since  ion-ion  collisions  at  high  densities  insure  that  indivi- 
dual nl  sublevels  will  equilibrate  to  their  statistically  weighted  densities  on  very  short  time  scales 
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(up  to  some  /ima„  that  is  where  A£'„„  becomes  small  and  A n - 1 transitions  become  as  likely; 
then  the  statistical  population  among  the  /-sublevels  is  re-distributed).  A study  was  also  made 
to  determine  whether  this  simplifying  approach  to  the  level  modeling  could  be  applied  to  heli- 
umlike ions.  It  was  found,  also,  that  ion-ion  collisions  and  opacity  effects  in  dense  plasmas  at 
high-temperature  can  successfully  re-distribute  the  populations  of  levels  having  principal  quan- 
tum number  greater  than  2,  such  that  the  / sub-levels  of  an  n level  are  maintained  at  their  sta- 
tistically weighted  densities  and  can  mathematically  be  represented  as  a single  degenerate  level. 
The  importance  of  the  environmental  effects  of  the  plasma  on  this  value  of  n.  however,  should 
be  re-emphasized,  and  any  plasma  environment  differing  significantly  from  the  one  considered 
here  must  be  re-assessed. 

A phenomenological  radiation  transport  scheme  has  been  implemented  in  our  plasma  radi- 
ation model  to  partially  account  for  the  effects  of  self-absorption  on  population  densities  and 
subsequent  radiation  emission.  An  attempt  to  compare  this  model  with  a spherical  ray-trace  cal- 
culation 24  is  underway  to  hopefully  provide  a means  of  assessing  the  capability  of  our  frequency 
diffusion  transport  as  well  as  to  improve  it. 

"Thin  plasma"  comparisons  between  this  model  and  a similar  rate  equation  scheme  have 
resulted  in  some  interesting  conclusions  concerning  the  completeness  of  an  atomic  structural 
model  for  a given  ion.  Good  agreement  was  obtained  between  the  two  calculations  for  popula- 
tion densities  of  the  heliumlike  ion  (where  level  structure  was  similar),  usually  within  15%. 
However,  for  the  hydrogenlike  ion.  only  estates  were  modeled  in  the  other  simulation,  whereas 
a complete  set  of  levels  up  to  n — 4 was  included  in  our  calculation.  In  addition,  no  collisional 
coupling  was  included  between  the  p states  of  the  other  model.  The  discrepancy  in  the 
predicted  values  for  populations  of  these  excited  levels  in  the  two  calculations  can  be  as  much 
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as  a factor  of  2.  Hence,  we  conclude  that  omission  of  the  neighboring  s,  d,  and  / angular 
momentum  levels  from  a structural  model  may  result  in  a significant  loss  of  accuracy  in  excited 
state  populations  due  to  collisional  and  radiative  coupling  to  these  missing  states.  Moreover, 
we  And  as  was  substantiated  by  Fig.  20,  that  in  failing  to  include  these  levels,  a serious  loss  of 
radiation  is  incurred  in  the  calculation  which  results  from  inter-level  radiative  transitions 
between  stales  greater  than  n - l and  for  which  Sn  * 0.  A/  - 1.  Our  simulation,  as  shown  in 
this  graph,  predicts  significant  radiation  from  these  transitions,  enhanced  by  the  reabsorption  of 
the  resonant  line  radiation. 

Comparisons  of  calculations  between  optically  thick  and  optically  thin  plasmas  have 
demonstrated  some  profound  effects  of  reabsorption  on  the  population  densities  (of  ground 
states  as  well  as  excited  states)  and  the  line  emission  from  optically  thirk  plasmas.  These  com- 
parisons have  shown  that  the  optical  pumping  of  excited  levels  due  to  trapped  radiation  and 
subsequent  strong  collisional  coupling  upward  can  cause  dramatic  increases  in  population  densi- 
ties of  higher  level  excited  states.  This  well-known  effect  of  opacity  on  the  radiation  can  result 
in  enhancement  of  thin  line  emission  often  by  as  much  as  an  order  of  magnitude  or  more,  as 
shown  in  Figs.  19  and  20  (see  specifically  the  We„.  and  He(2p-3d)  lines). 

In  addition,  we  have  demonstrated  that  the  increased  populations  of  the  excited  states  of 
ion  Z due  to  this  optical  pumping  can  cause  a similar  increase  in  the  ground  state  population  of 
the  Z + 1 ion  through  the  collisional  ionization  process.  This,  in  turn,  can  cause  increases  in 
populations  of  the  excited  states  of  the  Z + 1 ion  to  the  extent  that  this  increase  may  result  in 
such  a significant  enhancement  of  the  radiation  from  these  levels  that  it  will  overshadow  any 
decrease  (relative  to  the  thin  case)  one  would  expect  due  to  the  thickness  of  the  line.  The  star- 
tling result,  therefore,  is  that  an  optically  thick  line  may  emit  more  radiation  than  it  would  if  it 
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were  calculated  in  the  optically  thin  approximation,  as  does  the  H(L„ ) line  in  Fig.  19.  The 


neglect  of  this  possible  enhancement  could  have  a decidedly  detrimental  effect  on  the  diagnostic 


interpretation  of  spectral  line  emission  from  dense  laser-heated  plasmas  and  the  plasma  parame 


ters  which  are  inferred  from  their  ratios 
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Figure  2 — Argon  energy  level  diagram 
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Figure  8 — Electron  collision;*!  excitation  rale  coefficients 
for  heliumlike  aluminum  transitions 
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Entire  10  — Electron  collisinnal  excitation  rate  coefficients 
lor  hcliumliKe  alumimim  transitions 


Figure  13  - Electron  collisional  excitation  rate  coefficients 
for  heliumlike  aluminum  transitions. 
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